Q 2 evolution of parton distributions at small x 
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We investigate the Q 2 evolution of parton distributions at small x values, obtained in the case of flat initial 
conditions. The results are in excellent agreement with deep inelastic scattering experimental data from HERA. 



1. Introduction 

The measurements of the deep-inelastic scat- 
tering structure function F 2 in HERA Q have 
permitted the access to a very interesting kine- 
matical range for testing the theoretical ideas 
on the behavior of quarks and gluons carrying 
a very low fraction of momentum of the proton, 
the so-called small x region. In this limit one 
expects that non-perturbative effects may give 
essential contributions. However, the reasonable 
agreement between HERA data and the next-to- 
leading order (NLO) approximation of perturba- 
tive QCD that has been observed for Q 2 > lGeV 2 
(see the recent review in &) indicates that per- 
turbative QCD could describe the evolution of 
structure functions up to very low Q 2 values, tra- 
ditionally explained by soft processes. It is of 
fundamental importance to find out the kinemati- 
cal region where the well-established perturbative 
QCD formalism can be safely applied at small x. 

The standard program to study the small x 
behavior of quarks and gluons is carried out by 
comparison of data with the numerical solution 
of the DGLAP equations by fitting the param- 
eters of the x profile of partons at some initial 
Qq and the QCD energy scale A (see, for exam- 
ple, |j||l|])- However, if one is interested in ana- 
lyzing exclusively the small x region (x < 0.01), 
there is the alternative of doing a simpler analysis 
by using some of the existing analytical solutions 
of DGLAP in the small x limit (see [0 for re- 
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view). This was done so in Ref. §-§ where 
it was pointed out that the HERA small x data 
can be interpreted in terms of the so called dou- 
bled asymptotic scaling phenomenon related to 
the asymptotic behavior of the DGLAP evolution 
discovered in |gj many years ago. 

Here we illustrate results obtained recently in 
These results are the extension to the NLO 
QCD approximation of previous leading order 
(LO) studies |,§. The main ingredients are: 

1. Both, the gluon and quark singlet densities 
are presented in terms of two components ('+' 
and '— ') which are obtained from the analytical 
Q 2 dependent expressions of the corresponding 
('+' and ') parton distributions moments. 

2. The '— ' component is constant at small x, 
whereas the '+' component grows at Q 2 > Qq as 
~ exp (cr) , where 



a = 2y {d+s + D+p)lnx, 

and the LO term d + = — 12//3q and the NLO one 
D± = d± ± + d±[3 1 /[3 with d ±± = 412//(27A>). 
Here the coupling constant a s — a s /(47r), s = 
ln{a(Ql) / a(Q 2 )) andp = a{Ql) - a(Q 2 ), (3 and 
/3i are the first two coefficients of QCD /3-function 
and / is the number of active flavors. 

2. Basical formulae 

Thus, our purpose is to demonstrate the small 
x asymptotic form of parton distributions in the 
framework of the DGLAP equation starting at 
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some Q 2 with the flat function: 



fa(Q 



A„ 



(a = q,g), 



(1) 



where f a are the parton distributions multiplied 
by x and A a are unknown parameters that have 
to be determined from data. Through this work 
at small x we neglect the non-singlet quark com- 
ponent. 

In |^| an effective method to reproduce the x- 
dependence of parton distributions has been de- 
veloped. It is based on a separation of the sin- 
gular and regular parts of the exact solution for 
the moments of parton distributions and on the 
method to replace Mellin convolution by usual 
products The method allows in simplest way 
to reproduce the LO results jj| and to construct 
the ^-dependence of parton densities at NLO 0: 

fa(x,Q 2 ) = f+(x,Q 2 ) + f-(x,Q 2 ) and 
f~(x,Q 2 ) = A-(Q 2 ,Ql)exp(-d_(l) s 
-D_(l)p) + 0(x), 



f+(x,Q 2 ) = A+(Q 2 ,Q 2 )I (a)exp(-d + (l)s 
-D + (l)p) + O(p), 



(2) 



f+(x,Q 2 )=A+(Q 2 ,Q 2 ) 



(l-4(l)a(Q 2 )) 



xph(a) + 20a(Q 2 )I Q {a) 



xexp(-d + (l)s-D + {\)p) + O(p), 
F 2 (x,Q 2 ) =e(f q (x,Q 2 ) + ^fa(Q 2 )f g (x,Q 2 ) 



where are modified Bessel functions, which 

have z^-independent limit exp(a) at a — ► oo, 
e = (J2\ e i)/f is the average charge square for 
/ active quarks, p = a/2/ln(l/x) and the magni- 
tudes 

8 °MQ 2 ) 



At{Q 2 ,Ql) 



i 



81' 



A„ 



4r 
9 



1 



1 + 3(1 + -/)^) 



fa(Q 2 
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(3) 



81 J/ v u/ 81 
A-(Q 2 ,Q 2 ) = A g -A+(Q 2 ,Q 2 ), 

A+ = £ U g + , A- = A q - 20a(Q 2 )A+ 

s From now on, for a quantity k(n) we use the notation k 
for the coefficient in the front of the singular part when 
n — ► 1 and k(n) for the corresponding regular part. 



The regular parts of the terms d± and D± have 
the form n: 



d+(l) = l + — /, d_(l) 



d++(l) = - 36C3 + 33C2 
Po 
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where are Euler ^-functions. 
3. Slopes 

The behaviour of eqs. (0) can mimic a power 
law shape over a limited region of x, Q 2 : 

f a (x,Q 2 ) ~ ^-^(-.Q 2 ) and 
F 2 (x,Q 2 ) ~ x -*'/ 2 f {x,Q 3 ) 

The quark and gluon effective slopes \ e Jf are re- 
duced by the NLO terms that leads to the de- 
creasing of the gluon distribution at small x. For 
the quark case it is not the case, because the nor- 
malization factor A+ of the "+" component pro- 
duces an additional contribution undampening as 
-{Inx)- 1 ' 2 . 

The gluon effective slope \ e Jf is larger than the 
quark slope \ e q '' , that is in excellent agreement 
with a recent MRS and GRV analyses 
Indeed, the effective slopes have the asymptotical 
values (at large Q 2 ): 



\ e /f> as (x,Q 2 ) « p- 
\ e /f' as (x,Q 2 ) » 
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ihx(l/x) 
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\ e J 2 f > as (x,Q 2 ) « \ e /f> as (z,Q 2 ) + 
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where symbol « marks approximations ob- 
tained by expansions of modified Bessel func- 
tions I n (&). The slope A^/ 2 /,as (a;, Q 2 ) lies between 

4 The nonzero components of the singular parts were given 
in Introduction. 
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quark and gluon ones but closely to quark slope 
X e / f ' as (x,Q 2 ) (see also Fig. 2). 

Both slopes Xgff decrease with decreasing x. A 
x dependence of the slope should not appear for 
a parton density with a Regge type asymptotic 
(x~ x ) and precise measurement of the slope X e Jf 
may lead to the possibility to verify the type of 
small x asymptotics of parton distributions. 

4. Results of the fits 

With the help of the results obtained in the pre- 
vious section we have analyzed F^ HERA data at 
small x from the HI collaboration (first article in 
(D). In order to keep the analysis as simple as 
possible we have fixed ^y^( n f = 4) = 250 MeV 
which is a reasonable value extracted from the 
traditional (higher x) experiments. The initial 
scale of the parton densities was also fixed into 
the fits to Qq = 1 GeV 2 , although later it was 
released to study the sensitivity of the fit to the 
variation of this parameter. The analyzed data 
region was restricted to x < 0.01 to remain within 
the kinematical range where our results are accu- 
rate. Finally, the number of active flavors was 
fixed to /=4. 

Fig. 1 shows Fi calculated from the fit with Q 2 
> 1 GeV 2 in comparison with HI data. Only the 
lower Q 2 bins are shown. One can observe that 
the NLO result (dot-dashed line) lies closer to the 
data than the LO curve (dashed line). The lack 
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of agreement between data and lines observed at 
the lowest x and Q 2 bins suggests that the flat be- 
havior should occur at Q 2 lower than 1 GeV 2 . In 
order to study this point we have done the analy- 
sis considering Q 2 as a free parameter. Compar- 
ing the results of the fits (see ||) one can notice 
the better agreement with the experiment of the 
NLO curve at fitted Q 2 = 0.55GeV 2 (solid curve) 
is apparent at the lowest kinematical bins. 

Finally with the help of Eqs. (||) we have es- 
timated the F2 effective slope using the value of 
the parameters extracted from NLO fits to data. 
For HI data we found 0.05 < X e J 2 f < 0.30 - 0.37. 
The lower (upper) limit corresponds [] to Q 2 = 1.5 
GeV 2 (Q 2 = 400 GeV 2 ). The dispersion in some 
of the limits is due to the x dependence. Fig. 2 
shows that the three types of asymptotical slopes 
have similar values, which are in very good agree- 
ment with HI data (presented also in Fig. 2). 
The NLO values of A^ 2 h e between the quark 
and the gluon ones but closer to the quark slope 
X e q ff - as . These results are in excellent agreement 
with those obtained by others (see the review j|] 
and references therein). 



5. Conclusions 

We have shown that the results developed re- 
cently in Q have quite simple form and reproduce 



5 For smallO 2 we used the exact values of the slopes pre- 
sented in g. 
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many properties of parton distributions at small 
x, that have been known from global fits. 

We found very good agreement between our 
approach based on QCD at NLO approximation 
and HERA data, as it has been observed earlier 
with other approaches (see the review Q). Thus, 
the nonperturbative contributions as shadowing 
effects, higher twist effects and others seems to be 
quite small or seems to be canceled between them 
and/or with ln(l/x) terms containing by higher 
orders of perturbative theory (see discussion also 
in ]Io|). In our opinion, this very good agreement 
between approaches based on perturbative QCD 
and HERA data may be explained also by the fact 
that at low x values the real effective scale of cou- 
pling constant is like Q 2 /x c , where 1/2 < c < 1 
(see To clear up the correct contributions of 

nonperturbative dynamics and higher orders con- 
taining large ln(l/x) terms, it is necessary more 
precise data and further efforts in developing of 
theoretical approaches. 
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